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ABSTRACT 

To investigate the relationship between black holes and their host galaxies, many groups have used the width 
of the [O III] A5007 line as a substitute for the stellar velocity dispersion (cr*) of galaxy bulges. We directly test 
this assumption with a large and homogeneous sample of narrow-line active galactic nuclei from the Sloan Digital 
Sky Survey. We consider multiple transitions ([O II] A3727, [O III] A5007, and [S II] AA6716,6731) and various 
techniques for quantifying the line width in order to obtain a calibration between the gas velocity dispersion, dg, 
and cr*. We hnd that dg of the low-ionization lines traces as does dg for the core of [O III] after its asymmetric 
blue wing is properly removed, although in all cases the correlation between dg and d* has considerable scatter. 
While the gas kinematics of the narrow-line region of active galaxies are primarily governed by the gravitational 
potential of the stars, the accretion rate, as traced by the Eddington luminosity ratio, seems to play an important 
secondary role. Departures from virial motions correlate systematically with accretion rate. We discuss the 
implications of these results for previous studies that use [O III] line widths to infer stellar velocity dispersions in 
quasars and narrow-line Seyfert 1 galaxies. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: Seyfert 


1. INTRODUCTION 

The recent discovery of a tight correlation between stellar 
bulge velocity dispersion and central black hole (BH) mass (the 
ALbh - cr* relation) strongly indicates that BHs coevolve with 
their host galaxies (Gebhardt et al. 2000a; Eerrarese & Merritt 
2000; Tremaine et al. 2002). The growth and accretion histories 
of BHs appear to be inextricably linked with the mass assem¬ 
bly of galaxies. Thus, the evolution of BH mass density over 
cosmic time provides observational constraints on galaxy evo¬ 
lution. The Mbh - cr* relation allows us to track BH mass den¬ 
sity with time, provided we can measure a velocity dispersion 
for a given system. 

Unfortunately, the bulge stellar velocity disperion is often 
difficult to measure directly. In active galaxies the stellar ab¬ 
sorption features are quickly swamped by nonstellar emission 
from the nucleus. There is consequently a disconnect between 
the methods used to obtain BH masses for active galactic nuclei 
(AGNs) and quiescent galaxies. While locally it is possible to 
infer a BH mass function based on the distribution of stellar ve¬ 
locity dispersion (e.g., Yu & Tremaine 2002), at signihcant dis¬ 
tances we must rely on AGNs, which are significantly brighter 
than stellar light. The challenge remains to develop secondary 
estimators of BH masses in AGNs and assure ourselves that 
AGNs follow the same Mbh “ cr* relation. 

Luckily, there are a number of secondary techniques in devel¬ 
opment. Reverberation mapping provides a virial estimate of 
BH mass to within a factor that depends on the geometry of the 
broad-line region (BLR) (e.g., Kaspi et al. 2000; Peterson et al. 
2004). This factor can be empirically determined through com¬ 
parison of Mbh inferred from cr* and the reverberation mapping 
mass (Gebhardt et al. 2000b; Eerrarese et al. 2001; Nelson et al. 
2004; Onken et al. 2004). Reverberation mapping also provides 


an empirical calibration between AGN luminosity and BLR ra¬ 
dius, which may be combined with a BLR velocity to yield a 
virial mass estimate (e.g., Kaspi et al. 2000). This technique 
has been used widely in the literature (e.g., McLure & Dunlop 
2004) to study the cosmic evolution of BH mass density. 

In addition to secondary mass estimates for the BH, there is 
a secondary technique to estimate cr* using the velocity dis¬ 
persion of the ionized gas (cTg) from the narrow-line region 
(NLR). The NLR is in a privileged position in the galaxy; 
compact enough to be illuminated by the active nucleus and 
yet large enough to feel the gravitational forces of the bulge. 
Early studies of integrated line prohles determined that grav¬ 
ity, not nuclear activity, dominates the global kinematics of the 
NLR (Wilson & Heckman 1985; Whittle 1992a,b). Most re¬ 
cently, Nelson & Whittle (1996) hnd dg ~ cr* for a sample 
of 66 Seyfert galaxies. These studies all measure the width 
of [O III] A5007 because it is a strong and ubiquitous line. It 
seems that neither rotational support nor nuclear activity is im¬ 
portant compared to pressure support for the NLR, so that ccg 
can be used as a proxy for cr*. Of course, NLR gas velocity dis¬ 
persion, particularly of the [O III] line, is easier to measure for 
both active and distant objects. Many subsequent studies have 
used the full-width at half maximum (EWHM) of [O III] as a 
proxy for cr*. Nelson (2000) compares Mbh inferred from cTg to 
reverberation-mapped masses, while Boroson (2002) compares 
CTg to virial masses estimated using the line width-luminosity 
relation. Shields et al. (2003) compare the cosmic evolution 
of BH mass to the evolution of bulge velocity dispersion as in¬ 
ferred from CTg, while both Grupe & Mathur (2004) and Wang & 
Lu (2001) investigate the BH mass distribution of narrow-line 
Seyfert 1 galaxies (NLSls) using erg. 

Despite the widespread use of erg, caution is in order. While, 
to hrst order, the NLR kinematics seem to be determined by 
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gravity, the effects of the AGN are also apparent. Eirst of all, 
the [O III] line often has a broad blue wing (Heckman et al. 
1981; De Robertis & Osterbrock 1984; Whittle 1985a; Wilson 
& Heckman 1985 and references therein) and Whittle (1992b) 
finds that powerful, compact linear radio sources have anoma¬ 
lously broad [O III] lines. At least in some cases nuclear activity 
directly affects the NLR kinematics, and perhaps for these sys¬ 
tems a different emission line, one originating farther from the 
central source, would be a better tracer of cr*. It is also con¬ 
ceivable that rotation provides a significant source of support 
for the gas in some fraction of NLRs, thus lowering (jg relative 
to cr*. In short, it remains to be demonstrated with a large and 
homogenous data set that (jg truly traces cr*. The data for such a 
study are now readily available from the Sloan Digital Sky Sur¬ 
vey (SDSS; York et al. 2000; Stoughton et al. 2002). Using the 
database presented by Brinchmann et al. (2004a) we directly 
compare dg and cr* for a large sample of narrow-line AGNs 
from the SDSS. In addition to the conventional [O III] line, we 
consider the [Oil] AA3726, 3729 and [S II] AA6716, 6731 dou¬ 
blets. While these lines are problematic due to blending, and 
are typically weaker than the [O III] line, they are less contam¬ 
inated by strong asymmetries (as we show below). In addition, 
the [O II] line is observable to higher redshift. 

Ionized gas kinematics in the centers of galaxies is of broader 
interest than simply as a BH mass indicator. The structure and 
kinematics of the gas relative to the stars may provide clues 
about the dynamical history of the galaxy (e.g., Sofue & Ru¬ 
bin 2001). In many cases, it seems that the gas rotates more 
slowly than the stars, in both elliptical (e.g., Caldwell, Kirshner, 
& Richstone 1986; Cinzano & van der Marel 1994) and spiral 
galaxies (Eillmore, Boroson, & Dressier 1986; Kent 1988; Ko- 
rmendy & Westpfahl 1989; Bertola et al. 1995; Eisher 1997; 
Cinzano et al. 1999). As noted by many of these authors, this 
may indicate significant pressure support for the gas. On the 
other hand, Vega Beltran et al. (2001) find the gas in the cen¬ 
ters of ~ 20 disk galaxies is moving on circular orbits in a cold 
disk. Presumably, a variety of dynamical states are possible for 
the nuclear gas. It would be interesting to examine this ques¬ 
tion with a statistical sample such as SDSS. In particular, is 
there anything special about the dynamical state of the ionized 
gas in AGNs? 

We note that we are by no means the first to perform an anal¬ 
ysis of this type. Earlier studies have looked in detail at both the 
[O III] line shape as an indicator of the kinematics of the NLR 
and at the correlation of the line width with other nuclear and 
host properties (for a review, see Wilson & Heckman 1985 and 
references therein; Nelson & Whittle 1996). The new aspect 
of this work is the large, homogenous sample considered. It is 
also, to our knowledge, the first systematic investigation of the 
behavior of the [S II] and [O II] line profiles. Eurthermore, this 
is the first systematic investigation of NLR kinematics since the 
discovery of the Mbh - ct* relation. Not only does this provide 
new motivation for a better understanding of the cTg-cr* rela¬ 
tion, but it enables us to examine the kinematics of the NLR 
in relation to the physical properties of the AGN. In particular, 
we find that the accretion rate of material onto the BH partially 
determines the kinematics of the NLR. 

2. THE SAMPLE 


The Second Data Release (DR2) of the SDSS (Abazajian 
et al. 2004), with 375,000 galaxy spectra, enables us to com¬ 
pare nuclear gas and stellar kinematics with a large, homoge¬ 
neous sample. The reader is referred to the project website for 
more details, as well as a complete bibliography for the survey ^ 
Briefly, spectroscopic candidates are selected based on imaging 
in the u, g, r, i, and z bands (Eukugita et al. 1996; Smith et 
al. 2002; Strauss et al. 2002) with a drift-scan camera (Gunn 
et al. 1998). Spectra are acquired with a pair of fiber-fed spec¬ 
trographs. Each fiber subtends a diameter of 3", corresponding 
to ^5.4 kpc at z = 0.1. The spectra have an instrumental res¬ 
olution of A/AA w 1800 ((Jinst « 71 km s“'). Integration times 
are determined for a minimum signal-to-noise ratio (S/N) of 4 
at g = 20.2 mag. The spectroscopic pipeline performs basic im¬ 
age calibrations, as well as spectral extraction, sky subtraction, 
removal of atmospheric absorption bands, and wavelength and 
spectrophotometric calibration (Stoughton et al. 2002). 

Recently, the SDSS team at MPA/JHU made publicly avail¬ 
able their database of galaxy properties for 211,894 galaxies, 
including 33,589 narrow-line (Type 2) AGNs from DR2^. The 
sample selection and spectral treatment are discussed in detail 
in Kauffmann et al. (2003a). Like the main galaxy sample, 
the redshift distribution is strongly peaked at z « 0.1, and Ha 
is always present in the band. AGNs are selected based on 
their location in line intensity ratio diagnostic diagrams (Bald¬ 
win, Phillips, & Terlevich 1981; Veilleux & Osterbrock 1987; 
Kauffmann et al. 2003c) after the stellar continuum is mod¬ 
eled and removed using stellar population synthesis models. 
Their definition of AGNs is very inclusive. Only objects solidly 
on the stellar locus are excluded, while low-ionization nuclear 
emission-line regions (LINERs; Heckman 1980) and transition 
objects (Ho, Eilippenko, & Sargent 1993) are included. The 
database includes the line intensities used to make the selection. 
They also make available stellar velocity dispersions, measured 
using a direct-fitting algorithm that assumes Gaussian velocity 
profiles (Heckman et al. 2004). 

We selected our sample from this set of narrow-line AGNs. 
Using the tabulated emission-line ratios, we use the line-ratio 
diagnostics outlined in Ho, Eilippenko, & Sargent (1997a) to 
explicitly include only Type 2 Seyferts, excluding both LIN¬ 
ERs and transition objects. This guarantees that the emission 
lines are dominated by AGN light. As discussed in Ho (2004), 
the large physical scales subtended by the SDSS fibers very 
likely includes substantial host galaxy emission, which presents 
a source of confusion for interpreting LINERs and transition 
objects. We further require [O II], [O III], and [S II] to have a 
S/N per pixel greater than at least 3 in the line, and a measured 
stellar velocity dispersion. The mean equivalent widths (EWs) 
for the resulting sample are 18 for [O II], 20 for [O III], and 
4 for [S II], the weakest line. In §3.4 we limit our sample to 
objects having only the highest-EW lines to examine how the 
S/N affects our results. We exclude all objects with ct* below 
the nominal SDSS instrumental resolution of cTinst ~ 70 km s“' 
(as in Heckman et al. 2004). We further exclude a handful of 
objects with cr* > 400 km s“* as most likely unphysical (10 ob¬ 
jects had cr* w 600 km s“^)^ and 13 objects for which one of 
the line detections appeared marginal after manual inspection. 
This left a sample of 2050 objects. As described below, we also 
remove objects with double-peaked or unresolved line profiles. 


’ http://www.sdss.org/ 

^http://www.mpa-garching.mpg.de/SDSS/ 

^The distribution function of velocity dispersions for local early-type galaxies drops sharply for cr* > 400 km s“’ (Sheth et al. 2003). 
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leaving a final sample of 1749 objects. The median redshift of 
the sample is z « 0.1, the same as the SDSS in general. 

Throughout we assume the following cosmological parame¬ 
ters to derive distances: Ho = lOO/i = 72 km s“' Mpc“', = 
0.3, and ffA=0.7. 

3. DATA ANALYSIS 

We wish to compare the stellar and gas line widths for the 
entire AGN sample. Since the stellar velocity dispersions have 
been measured, we need only remove the stellar continuum and 
then measure the gas line widths. Below we describe each of 
those steps. 

3.1. Continuum Subtraction 

In order to measure reliable line widths, the underlying stel¬ 
lar continuum must be removed. Continuum variations can 
hide or mimic broad line wings and skew relative line inten¬ 
sity measurements. Since the MPA/JHU database does not pro¬ 
vide continuum-subtracted spectra, we have to repeat this step 
of the analysis. Fortunately, the SDSS survey provides a large 
and homogeneous library of pure absorption-line galaxies for 
use in continuum subtraction. We use a code that employs prin¬ 
cipal component analysis, kindly provided by Lei Hao (Hao et 
al. 2005). A library of hundreds of spectra is projected onto 
an orthogonal basis of eigenspectra. Typically eight eigen- 
components are sufficient to provide an adequate model of the 
stellar spectrum. A power-law continuum and an A star spec¬ 
trum may be included if necessary. 

3.2. Profile Fitting 

There are many possible ways to parameterize emission-line 
profiles. Some of the most common include the FWHM and the 
second moment of the line, defined to be 

a^ = {c/X,f j {\-Kffxd\/ J fxdX. 

The line centroid is Ac, f\ is the flux density in the continuum- 
subtracted spectrum, and a is in units of km s“'. For a pure 
Gaussian the second moment is equal to a and FWHM = 
aj2.35. (The second moment is directly comparable to the stel¬ 
lar velocity dispersion measurements, which assume a Gaus¬ 
sian broadening function.) If to this original Gaussian we added 
a broad, low-contrast Gaussian to represent the line wings, the 
FWHM would increase slightly, by an amount depending on 
the contrast between the core and the wings. The second mo¬ 
ment, on the other hand, is more sensitive to the presence of 
wings, and becomes significantly broader. For this reason, we 
get more information on the full line profile from the second 
moment. Unfortunately, measurements of the second moment 
can be very sensitive to the S/N of the spectrum. Therefore, 
we first fit a multi-component Gaussian to the line profile, and 
then measure the second moment from the model profile, as 
described below. This technique has the further benefit of pro¬ 
viding a decomposition of the core and wing components, so 
that we can examine the distribution of core and wing widths 
independently. 

After much experimentation, we find that the most reliable 
method, applicable to the majority of objects, is to fit each 


profile with two Gaussians, one representing the line core and 
the other representing the line wing. The second component 
is fixed initially to be blueward of the central component, al¬ 
though in a subset of cases (~ 6%), a redshifted component was 
required for an acceptable fit. If the fit with a second component 
improves significantly (at least 20% increase in x^) over the 
single-component fit, it is kept. We use this model to simulta¬ 
neously fit the [O III] AA4959, 5007 lines, with their separation 
fixed and their relative strength constrained to the theoretical 
ratio of 3:1. Likewise, the [S II] AA6716, 6731 lines are fixed to 
the same profile"*, with their separation fixed to the laboratory 
value; the ratio of the line intensities varies with density, and 
is allowed to vary in the fit. The [O II] AA3726, 3729 lines be¬ 
have similarly, but the decomposition is more difficult because 
the expected line widths are comparable to the pair separation 
(2.8 A or 225 km s“'). For this reason, we allow the line ratio 
to vary between 0.2 and 1.4, the expected range for the typical 
densities observed in the NLR (e.g., Osterbrock 1989). Sample 
fits are shown in Figure 1. 

As in previous studies, the majority of the line profiles are 
smooth, but a small percentage are lumpy, and some have two 
distinct peaks with comparable flux and line width. We find 31 
of the 2050 (~ 1%) objects with clearly defined double peaks 
(Fig. 1). While we show below that in general the NLR is not 
rotationally dominated, we suspect that in these cases the NLR 
gas is predominantly organized into a disk. Of course, bipolar 
outflows may also be responsible for the double-peaked lines. 
In any case, these objects are such clear outliers that we exclude 
them from further consideration. 

3.3. Instrumental Resolution 

The observed a is the convolution of the true line width and 
the instrumental response. To first order, the intrinsic line width 

can be approximated by cr = *he case of the 

SDSS data, (Tinst varies not only as a function of wavelength, 
but also the location of the object on the plate and the tempera¬ 
ture on the night of the observations. For this reason, the SDSS 
pipeline measures the instrumental response using an arc lamp 
spectrum, and returns the resolution at every pixel. The mean 
values for this sample are 74 ± 8 km s“' for [O II], 60 ± 5 km 
s“' for [O III], and 56 ± 5 km s“' for [SII]. When we remove the 
effects of instrumental resolution, some of the objects become 
unresolved. Measurements below the resolution limit are not 
reliable, and so we remove all objects for which one or more 
of the lines in question is unresolved, leaving a total of 1749 
objects. 

3.4. Errors and Uncertainties 

The average formal error in (jg is ~ 5%. However, this esti¬ 
mate neglects systematic effects, the two most serious of which 
are uncertainties in the continuum subtraction and component 
decomposition. The latter is only a concern in the case of 
[O III], the only line that often requires a two-component fit. 
In some cases, there is ambiguity as to the exact decomposition 
between these two. However, we do not believe errors in de¬ 
composition introduce any systematic bias. There are publicly 
available error estimates^ for the stellar velocity dispersions we 


■^If the NLR is stratified in density, the width of each component of the doublet can be different, since their critical densities are different. But the profile differences 
are expected to be small and hard to discern with data of moderate S/N. 

^http;//spectro.princeton.edu/ 
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i^IG. 1.— Sample fits for blue-asymmetric (top), red-asymmetric {middle), and double-peaked {bottom) objects. The ordinate of the plots are in units of 10' 
^ cm“^ The ordinate scales are 1:1:0.33 {top), 1:3:1 {middle), and 1:1:0.2 {bottom). Note the hint of double-peakedness in the [S ii] profile {bottom). 
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Table 1. Statistical Results 


Sample 

(1) 

Number 

(2) 

(3) 

(4) 

('^[OIII]/o'*) 

(5) 

(6) 

/o-*) 

All 

1749 

1.13 

± 

0.38 

1.11 

± 

0.35 

1.34 

± 

0.66 

1.00 

± 

0.35 

High EW 

657 

1.19 

± 

0.45 

1.15 

± 

0.38 

1.40 

± 

0.53 

0.98 

± 

0.35 

Low 2 : 

736 

1.17 

± 

0.39 

1.16 

± 

0.37 

1.43 

± 

0.75 

1.01 

± 

0.35 

Moderate 2 : 

957 

1.11 

± 

0.37 

1.09 

± 

0.34 

1.28 

± 

0.57 

1.00 

± 

0.35 

High 2 : 

957 

0.99 

± 

0.35 

1.00 

± 

0.36 

1.19 

± 

0.53 

0.95 

± 

0.27 

Low Ha 

872 

1.13 

± 

0.38 

1.13 

± 

0.38 

1.28 

± 

0.67 

1.04 

± 

0.35 

High Ha 

877 

1.14± 

0.38 

1.10 

± 

0.32 

1.40 

± 

0.63 

0.97 

± 

0.34 

No radio 

1475 

1.12 

± 

0.37 

1.10 

± 

0.34 

1.29 

± 

0.59 

1.00 

± 

0.34 

Low radio 

119 

1.27 

± 

0.43 

1.34 

± 

0.47 

1.80 

± 

0.99 

1.06 

± 

0.41 

High radio 

133 

1.20 

± 

0.43 

1.12 

± 

0.34 

1.56 

± 

0.76 

1.01 

± 

0.41 

Lobes 

23 

0.99 

± 

0.40 

0.93 

± 

0.32 

1.00 

± 

0.44 

0.75 

± 

0.22 


Note. — Samples High EW: EW([Oii])> 25 A, EW([Oiii]) > 20 A, and EW([Sii]) 
> 5 A, data from Brinchmann et al. 2004a. Low z: z < 0.1. Moderate z: 0.1 < z < 0.2. 
High z: z > 0.2. Low Ho: Lhq < erg High Ho: Lhq > erg 8“^. 

No radio: no detection in FIRST. Low radio: P 20 cm < 10^^'® W Hz“^. High radio: 
P 20 cm > 10^^ ® W Hz“^. Lobes: Spatially extended in FIRST. 


are using. The median error is ~ 8%, or ~ 12 km s“' at the me¬ 
dian cr* w 144 km s“'. These uncertainties are small compared 
to the scatter in (crg/cr*), the quantity central to our analysis in 
this paper. 

We can investigate the extent to which the scatter is due to 
measurement error by limiting the sample to a subset of high- 
EW objects for which the formal measurement errors are small. 
We limit the sample to the strongest ^ 500 objects for each 
line. We find no significant difference in either the mean value 
or the scatter around (jg (see entries in Table 1). This suggests 
that the scatter is dominated by intrinsic variations, rather than 
measurement errors. In the Appendix we specifically examine 
the scatter introduced by rotational broadening. We find that it 
is also insufficient to account for the observed scatter. 

4. RESULTS 

4.1. Primary Driver: Gravitational Potential of the Bulge 

Our findings are summarized in Table 1, which shows 
(cTg/cr*) for each transition. Our primary result is that the line 
widths of [S II] and [OII] track the stellar velocity dispersion in 
the mean, albeit with substantial scatter. The width of [O III], 
on the other hand, is on average significantly broader than ct*; 
Pioiii] cannot be used directly as a tracer of cr*. However, as we 
discuss below, the core of [O III] does track ct*. 

The different behavior of the low-ionization lines compared 
to [O III] can be understood in terms of well-known variations 
in line shape as a function of critical density and ionization po¬ 
tential within the NLR (e.g., Pelat, Alloin, & Fosbury 1981; Fil- 
ippenko & Halpern 1984; Filippenko 1985; De Robertis & Os- 
terbrock 1986). The NFR is stratified in density such that lines 
of higher critical density originate from closer to the nucleus. 
Since the critical density of [O III] A5007 is higher than those 
of [Oil] A3727 and [S II] AA6716,6731 (Osterbrock 1989), 
in a density-stratified NFR we expect [O III] to be broader 


than [O II] and [S II], whereas [O II] and [S II] should share 
similar line widths. As expected, we find that (cr[oui]/o'[oii]) 
and (cr[oui]/cr[su]) « 1.2 ±0.5, while (<T[sii]/cr[oii]) « 1.0 ±0.3. 
This suggests that the [O III] emission originates from closer to 
the nucleus than the other lines. 

As noted, our results suggest that the kinematics of low- 
ionization gas in the NFR is dominated by gravity, such that the 
velocity dispersion in the gas and stars are comparable. On the 
other hand, the [O III] line width is significantly broader than 
(j*, suggesting that the active nucleus plays a significant role in 
governing the kinematics of the [O III]-emitting gas. Taken at 
face value, this finding seems to contradict the conclusions of 
Nelson & Whittle (1996), who have argued that the [O III] line 
width is largely dominated by gravity. Since the [O III] profile is 
invariably asymmetric, characterized by a central core compo¬ 
nent plus a high-velocity blue wing, it is reasonable to suppose 
that the line core should have the closest affinity to any grav¬ 
itationally dominated component. Using our two-component 
Gaussian fits to the [O III] profiles, we can directly compare the 
core velocity dispersion with cr*. The result of this test is shown 
in the last column of Table 1. Indeed, we find that when the blue 
wing is removed, the core of the [O III] line statistically traces 
again with significant scatter. This suggests that the low- 
velocity gas in the line core is gravitationally dominated, while 
the wing is more strongly influenced by the active nucleus. 

4.2. Secondary Drivers 

We now look for secondary parameters to explain the de¬ 
viations of (Tg for the [O III] line compared to ct*, which are 
plausibly related to either the active nucleus or the environ¬ 
ment. Fooking at the correlation of these deviations with vari¬ 
ous indicators of nuclear activity may reveal useful information 
about the interaction of the nucleus with the NFR. We con¬ 
sider host galaxy morphology, local environment, star forma- 
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tion rate, and bulge stellar velocity dispersion as global vari¬ 
ables. In terms of nuclear effects, we look at radio power, 
AGN luminosity, and Eddington ratio [Lboi/EEdd, where LEdd = 
1.26 X 10^* erg s“^ (Mbh/Mq)]. In particular, we will exam¬ 
ine the correlation of each of the above properties with Act 
= log CT[oiii] - log CT*. In Figure 2 we show Act plotted against 
each variable, with binned values overplotted to highlight any 
underlying trend. The large, binned points represent the mean 
and standard deviation in the mean plotted at the bin center. 

Of additional interest is whether the characteristics of the 
blue wing correlate with any of the properties listed above. In 
order to investigate this question we looked at three wing diag¬ 
nostics: the width of the blue wing, the flux ratio between the 
core and the wing, and the wavelength shift between the core 
and the wing. We found no strong correlations between any of 
the above quantities with either the nuclear or global variables. 

4.2.1. Global Properties 

Both gas and stellar kinematics in the central regions of 
galaxies may depend on Hubble type. We cannot derive robust 
morphologies directly from the SDSS data because the reso¬ 
lution and depth are insufficient. However, Shimasaku et al. 
(2001) and Strateva et al. (2001) find that the color, concentra¬ 
tion index, and spectral decomposition provided by the SDSS 
pipeline all correlate well with standard morphological classifi¬ 
cations. We use two different integrated properties to estimate 
the morphological distribution of the sample. The “eClass” is 
a spectroscopic classification derived from the principal com¬ 
ponent analysis decomposition of each galaxy spectrum, and 
ranges from -0.35 for early-type galaxies to h- 0.5 for late-type 
galaxies. The “concentration index,” the ratio of the radii en¬ 
closing 90% and 50% of the light in the r band, describes the 
compactness of the stellar emission. Concentrations indices 
larger than 2.6 correspond to early-type galaxies. The mean 
concentration index of 2.65 and mean eClass of-0.05 are both 
consistent with a population of early-type galaxies (Kauffmann 
et al. 2003c). We see no trend with Act in either of these pa¬ 
rameters (see Fig. 2a for Act vs. eClass). 

Nelson & Whittle (1996) find evidence for broadened line 
widths in noticeably disturbed systems. While we cannot use 
the SDSS imaging to reliably gauge disturbance levels, we can 
look for systematic trends with local galaxy density. Kauff¬ 
mann et al. (2004) find an increase of nuclear activity and star 
formation for AGNs in overdense regions. We might, there¬ 
fore, expect some correlation between environmental density 
and line width. However, we find no evidence for a trend. The 
luminosity-weighted neighbor counts are plotted against Act 
for those objects with measured neighbor counts in Figure 2b. 

Depending on the origin of the ionized gas, star formation 
(particularly in the nuclear regions) could affect the gas kine¬ 
matics. Recent studies, such as that of Kauffmann et al. (2004), 
suggest a link between star formation and nuclear activity. If 
this is the case, then stellar winds from massive stars or super¬ 
nova feedback may impact the dynamical state of the nuclear 
gas. To investigate this possibility, we use the total star forma¬ 
tion rate per unit stellar mass (SFR/M*) from Brinchmann et al. 
(2004b), which is derived for the AGN sample using a relation 
between the 4000 A-break and SFR in a sample of star-forming 
galaxies. In Figure 2c we plot log SFR/M* against Act; there is 


no significant trend apparent. 

By far the strongest trend is between Act and ct* (Fig. 2d), 
with a Spearman rank correlation coefficient of p = -0.32, and a 
probability of Pnuii < 10““* for the null hypothesis of no correla¬ 
tion. Massive galaxies with large values of ct* on average have 
Act « 0, whereas less massive galaxies with smaller values of 
CT* systematically have Act > 0. This trend with ct* directly 
translates into a trend with BH mass, via the Mbh - cr* relation. 

4.2.2. Radio Power 

Many lines of evidence suggest that radio jets influence the 
physical and dynamical state of the NFR. Wilson & Willis 
(1980) found a correlation between [O III] line width and ra¬ 
dio power. There is also a well-known relation between [O III] 
luminosity and radio power (de Bruyn & Wilson 1978; Ho & 
Peng 2001). Whittle (1985b,1992a,b) finds evidence that strong 
(Ti 415 MHz ^ 10^^'^ W Hz“'), “linear”® radio sources have sys¬ 
tematically broader [O III] lines. The same effect is seen by 
Nelson & Whittle (1996), but the trend is weaker. It seems 
that the jet plasma interacts with and accelerates the NFR, thus 
boosting the line width. We stress that this finding holds only 
for compact linear radio sources. In contrast. Smith, Heckman, 
& Illingworth (1990) do not find supervirial [O III] line widths 
among powerful radio galaxies with kpc-scale jets. 

In light of the above results, we extracted radio fluxes from 
the Faint Images of the Radio Sky at Twenty-cm (FIRST; 
Becker, White, & Helfand 1995) survey, which is largely co¬ 
incident with the SDSS footprint. FIRST is a 20 cm survey 
with a sensitivity of ~ 1 mJy and a resolution ~ 5". All but 
279 of the sources in our sample have FIRST measurements. 
The SDSS pipeline only considers matches within 2" in order 
to minimize confusion (Ivezic et al. 2002). However, a small 
fraction of the sources will have extended sources associated 
with them as well. We therefore search the FIRST database 
with a 60" radius for extended sources, and find 25. The radio 
luminosities are K-corrected to 20 cm assuming /v fx 
which is the median spectral slope found by Ho & Ulvestad 
(2001) for nearby Seyfert galaxies. Unfortunately, while the 
FIRST database is homogeneous and complete, the resolution 
is not sufficient to resolve nuclear morphology. Nevertheless, 
for completeness we consider possible trends with Act. 

Our results are ambiguous (see Fig. 2e). As can be seen from 
Table 1, the presence of core radio emission seems to have no 
impact on the observed line widths. On the other hand, ex¬ 
tended radio sources appear to have marginally narrower line 
widths. We attribute the ambiguous results to the low sensitiv¬ 
ity and resolution in the radio data. The correlation previously 
discovered holds for linear radio cores only. With 5" resolution, 
we cannot discriminate between different core morphologies. 
Furthermore, the sensitivity limit of 1 mJy corresponds to 
PiQcm = 2.5 X 10^^ W Hz“', at the “powerful” limit considered 
by Whittle (1992b). 

4.2.3. AGN Luminosity and Eddington Ratio 

Perhaps the most basic question to ask is whether the energy 
output of the active nucleus is directly related to Act. We use 
emission-line (Ha and [O III]) luminosity as a surrogate for the 
AGN luminosity (Adams & Weedman 1975; Yee & Oke 1978; 
Yee 1980). We derive Lhq and L[om] from the fluxes given in 


®The “linear” classification for the above studies comes from the high-resolution (072-076 or 80^00 pc) study of Ulvestad & Wilson (1984). Linear sources 
have extended radio emission on sub-galactic scales (tens to hundreds of parsecs), as opposed to sources classified as “diffuse”, “slightly resolved”, or “unresolved.” 
The linear morphology suggests the presence of compact jets, not to be confused with the extended jets and radio lobes found in radio galaxies. 
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the database of Brinchmann et al. (2004a), corrected for Galac¬ 
tic extinction using the dust maps of Schlegel, Einkbeiner, & 
Davis (1998). We find no evidence for a strong trend with ei¬ 
ther quantity (see Eig. 2/for L[oiii]), at least for the relatively 
limited range of luminosities probed here. This result does not 
change for a high-EW subsample. 

In addition to the luminosity of the AGN, we can ask whether 
the Eddington ratio affects the NLR kinematics. The Edding¬ 
ton ratio, proportional to the accretion rate per unit (BH) mass, 
characterizes the extent to which radiation pressure competes 
with gravity in the nucleus. Therefore, we might expect objects 
with higher Eddington ratios to show larger nonvirial motions. 

We do not have direct measurements of Eddington ratios for 
our sample, but we can estimate them by using line luminosi¬ 
ties as a substitute for Lboi and ct* as a surrogate for BH mass 
(and hence LEdd)- We stress, however, that the estimated val¬ 
ues of Lboi/TEdd are highly uncertain. The conversion from the 
observed luminosity in a specific band (e.g., optical) to bolo- 
metric luminosity introduces the largest uncertainty. Eor Type 2 
objects, the problem is compounded by the fact that the nonstel- 
lar (AGN) continuum is not directly observable in the optical. 
Instead, we will make use of the extinction-corrected [O III] lu¬ 
minosity to estimate Lboi, using the bolometric correction for 
[O III] luminosity from Heckman et al. (2004). This is derived 
in two steps. Eirst, the [O III] luminosity is correlated with op¬ 
tical AGN continuum luminosity for a sample of Type 1 ob¬ 
jects in SDSS (Kauffmann et al. 2003c; Zakamska et al. 2003). 
Then a standard bolometric correction from optical continuum 
to bolometric luminosity (from Marconi et al. 2004) is em¬ 
ployed. The final derived relation becomes Lboi ~ 3500L[oiii], 
with an uncertainty of 0.4 dex. We note that this estimate agrees 
well with the [O III] bolometric correction derived by L. C. Ho 
(2005, in preparation), established in a completely independent 
way using broad-band spectral energy distributions. The BH 
mass is calculated using the measured cr* and the Mbh - cr* re¬ 
lation of Tremaine et al. (2002), which formally has a scatter of 
less than 0.21 dex. 

Notwithstanding the large uncertainty in bolometric correc¬ 
tion, we find a highly significant correlation between Act and 
Tboi/TEdd (Fig. 3a). The [O III] line width is most supervirial 
in objects at the highest Eddington ratio. Formally, the Spear¬ 
man rank correlation coefficient is p = 0.46, with a probability 
of Lnuii < 10““* for rejecting the null hypothesis of no correla¬ 
tion. For the purpose of comparing with previous work, we 
also calculate Act using FWHM/2.35 rather than CTg. Although 
we have argued that the FWHM is an inferior measure of line 
width, there are occasions in which only FWHM measurements 
are available in the literature. The result is shown in Figure 
3L. While the slope is shallower, the trend persists (p = 0.36, 

Pnull < 10-“*). 

We have seen that Act is also significantly correlated with ct*, 
or, equivalently, BH mass (§4.2.1). Since Lboi/LEdd depends on 
BH mass, we must consider whether BH mass or Lboi/LEdd is 
really the fundamental parameter driving Act. Because the cor¬ 
relation of Act with Lboi/LEdd (p = 0.46) is stronger than that 
with BH mass (p = -0.32), Lboi/LEdd appears to influence Act 
more strongly than BH mass. In principle, the strongest corre¬ 
lation could come from some combination of the BH mass and 
Lboi/LEdd- However, we find that including the BH mass in a 
partial correlation analysis (Akritas & Siebert 1996) does not 
improve the correlation. We therefore conclude that Lboi/LEdd 
is the primary physical parameter driving Act. 


The solid line in Figure 3a represents our fit to the data. Us¬ 
ing an ordinary least-squares regression with log Lboi/LEdd as 
the independent variable, we find 

Act = (0.131 ± 0.006) log Lboi/LEdd+ (0.202 ± 0.006). (1) 

To better highlight the trend, we also binned the data in bins 
of A(log Lboi/LEdd) = 0.5 dex, which yields the following fit 
(dashed line): 

Act = (0.113 ± 0.005) log Lboi/f-Edd + (0.192 ± 0.008). (2) 

If we calculate Act using CTg = FWHM/2.35 (Fig. 3b), then the 
corresponding relations for the unbinned and binned data are, 
respectively. 

Act = (0.072 ± 0.005) log Lboi/LEdd + (0.080 ± 0.005) (3) 

and 

Act = (0.065 ± 0.002) log Lboi/LEdd+ (0.077 ± 0.004). (4) 

It is interesting to examine Act at the extreme Eddington ra¬ 
tios. For objects at the lowest Eddington ratios (log Lboi/LEdd 
^-2 for our sample) we find line widths that are subvirial 
(Act < 0). This finding suggests that quiescent (inactive) 
early-type galaxies (most closely approximated by the lowest- 
Eddington ratio systems) have subvirial gas velocities. In sup¬ 
port of this finding, Vega Beltran et al. (2001) find the CTg < ct* 
in early-type spirals, where the gas is at least partially rotation- 
ally supported (in contrast with late-type spirals, where both 
gas and stars are rotationally supported). In light of this issue, 
it is interesting to note that the lobe-dominated radio sources, 
which have been suggested to be preferentially low-Eddington 
ratio objects (e.g., Rees et al. 1982), also appear to have sub¬ 
virial gas velocities (§ 4.2.2). We do not find a significantly 
low distribution of Eddington ratios among the extended radio 
sources in our sample, but the statistics are poor. It would be 
interesting to extend our study to a larger sample of extended 
radio sources. 

At the other extreme, for log Lboi/LEdd ^ 1, our results im¬ 
ply that CTg systematically overestimates ct* by ~ 0.2-0.3 dex. 
This strongly suggests that studies that use CTg as a proxy for ct* , 
such as that by Shields et al. (2003) or Grupe & Mathur (2004), 
may systematically overestimate ct*. In both cases we expect 
the objects to be preferentially accreting at a high rate. Shields 
et al. (2003) concentrated on high-luminosity sources (QSOs) 
at high redshift, which would tend to be biased toward higher 
Eddington ratios. The study of Grupe & Mathur (2004) is con¬ 
cerned with NLS1 s, which are thought to be preferentially high- 
Eddington ratio objects (e.g.. Pounds, Done, & Osborne 1995). 
Grupe & Mathur (2004; see also Mathur, Kuraszkiewicz, & Cz¬ 
erny 2001), using FWHM[oiii] as a surrogate for ct*, claim that 
NLS Is have BHs that are undermassive relative to their bulge 
masses. Most or all of the trend discussed by these authors may 
be attributed to the systematic broadening of [O III] at the high¬ 
est Eddington ratios. 

In Figure 4 we show a modified version of Figure 1 from 
Grupe & Mathur (2004). We have reproduced the figure ex¬ 
actly in open symbols, and then overplotted the same data with 
the [O III] widths corrected using our derived calibration (Eq. 
3) in filled symbols. Again, while we do not advocate the 
use of FWHM, in this case it is the only available measure 
of line width. The squares are “normal” broad-line Seyfert 
Is, while the circles are NLS Is, selected using the Grupe & 
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Fig. 3.— (a) Act = log(j[oiii] -logcr* plotted against the Eddington ratio, Eboi/^Edd- Bolometric luminosity is estimated using the bolometric correction from 
Heckman et al. (2004); see text. Points are the SDSS sample. Large points are the mean values plotted in the center of the bin, and error bars are the standard 
deviation in the mean. Solid line shows the fit to data in Equation 1; dashed line shows the fit to the binned data in Equation 2. (b) Same as (a), but using cTg = 
FWHM/2.35. Boxes are the Type 1 objects from the Nelson & Whittle (1996) sample. The solid and dashed lines represent Equations 3 and 4, respectively. 


Mathur (2004) criterion of FWHM(H/?) < 2000 km The 
Mbh values are derived from FWHM(H/3) and L 5100 as tabu¬ 
lated in Grupe et al. (2004; corrected to our cosmology), using 
the line width-luminosity-mass scaling of Kaspi et al. (2000). 
The overplotted line represents the Tremaine et al. (2002) pa¬ 
rameterization of the Mbh-ct* relation. The Tboi/TEdd values 
used for the correction are derived using the Lboi values pro¬ 
vided in Grupe et al. (2004), estimated from their optical and 
soft X-ray data. The mean Eddington ratio for the broad-line 
Seyfert Is in this sample is log Lboi/Tedd = -0.74, while for 
the NLSls it is log Lboi/Tedd = 0.23, in support of the idea 
that NLS1 s tend to be high-Eddington ratio objects. This trans¬ 
lates into (Act) = 0.03 for the broad-line objects, as compared 
to (Act) =0.1 for the NLSls. As shown in Eigure 4, the ar¬ 
gument that the NLSls lie systematically below the Mbh-ct* 
relation is significantly weakened. To test this statistically, we 
compute the perpendicular distance (d±y of each point from 
the Tremaine et al. (2002) line. After our correction, the two 
populations are no longer statistically different in d± accord¬ 
ing to a Kolmogorov-Smirnov (K-S) test (Fks = 25%). In other 
words, the case that the BHs of NLSls are undermassive with 
respect to their bulges is no longer statistically supported. 

We are wary to overinterpret Eigure 4, since it is based on 
our statistical correction for Act and on a bolometric correc¬ 
tion, both of which inherently have a large scatter. In actuality, 
the best way to establish the location of NLSls relative to the 
Mbh - cr* relation is to obtain direct measurements of ct* for 
them. Barth, Greene, & Ho (2005), in a systematic investiga¬ 


tion of CT* for the NLSls selected from the sample of Greene & 
Ho (2004), show that these objects follow theMBH-cr* relation. 


5. DISCUSSION 
5.1. Seyfert Type 

We must examine the extent to which the results presented 
here are applicable to broad-line (Type 1) AGNs, since it is 
precisely in these objects that narrow emission-line widths are 
used to estimate ct*. According to the unification hypothesis 
of Seyfert galaxies, the differences between Type 1 and Type 2 
sources arise from orientation alone (e.g., Antonucci & Miller 
1985), such that the BLR in Type 2 AGNs is hidden from di¬ 
rect view by an obscuring torus. In this picture, the NLR is 
an extended, roughly isotropic region extending far beyond the 
torus. Thus, our results should apply to both AGN types. On 
the other hand, it is possible that the fractional incidence of blue 
wings depends on the obscuring geometry, which in turn would 
depend on Seyfert type. This could result, for instance, if the 
torus is partially responsible for collimating the ionizing radia¬ 
tion field (see, e.g., Tadhunter & Tsvetanov 1989). That said, to 
first order at least, our results should apply to Type 1 and Type 
2 objects. We offer two lines of evidence to support this. 

Eirst, if we consider well-studied, unbiased samples of 
Seyfert galaxies, such as that derived from the Palomar survey 


'If dx is the difference between the observed [Oiii] line width and the Tremaine et al. (2002) prediction based on the Mbh measurement, and dy is the difference 
between the observed Mbh ^nd that predicted by the measured [Oiii] line width, then d± = dx dy I yjd^ +d^. We (arbitrarily) define d± > 0 if the point lies below the 
Tremaine et al. (2002) line. We choose this statistic because it does not assume an independent variable. 
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Fig. 4. — Mbh vs. [O m] line width for the Grupe & Mathur (2004) sam¬ 
ple. We use the FWHM in this case because it is the quantity published. Open 
symbols are a reproduction of Figure 1 in Grupe & Mathur (2004), while filled 
symbols show the data for c7[oin] corrected for the Eddington ratio according 
to Equation 3. Circles are NLSls; boxes are broad-line Seyfert Is. All mea¬ 
sured quantities were taken from Grupe et al. (2004). The Mbh values come 
from the line width-luminosity-mass scaling relation of Kaspi et al. (2000). 
The solid line represents the Mbh “ cr* relation from Tremaine et al. (2002). 

(Ho et al.l997a,b), we see that the distribution of bulge proper¬ 
ties and nuclear emission-line widths are comparable for both 
Seyfert types. As shown in Table IB of Ho et al. (2003), there 
is no statistical difference between the Type 1 and Type 2 ob¬ 
jects in terms of Hubble type, absolute magnitude, and bulge 
luminosity, implying identical distributions of cr* between the 
two groups. At the same time, the distribution of NLR line 
widths (see Table 2B of Ho et al. 2003) are also statistically 
equivalent. The similarities in both bulge and nuclear gas prop¬ 
erties strongly suggest that Type 1 objects will show the same 
distribution in (dg/cr*), as that seen here for Type 2 sources. 

Second, we can directly test the hypothesis that (cTg/cr*) is 
invariant with respect to Seyfert type by using the study of Nel¬ 
son & Whittle (1996), which provides cr* and dg for 22 Type 1 
and 44 Type 2 objects (these are the 66 objects with a tabulated 
[O III] line width and stellar velocity dispersion). Again, only 
FWHM measurements are available. We compare FWHM[oiii] 
with d* for these objects in Figure 5. The two samples largely 
appear to overlap. For Type 1 objects, (dg/d*)= 1.00 ±0.46 
while for Type 2 objects (dg/d*)= 1.22 ±0.78. We use the 
K-S test to compare the distribution of dg/d* for each group 
and find that there is no statistically significant difference be¬ 
tween the two; the probability of rejecting the null hypothesis 
that the two groups are drawn from the same parent popula¬ 
tion is Fks = 13%. When we remove the four obvious outliers, 
which are the powerful linear radio sources discussed by Whit¬ 
tle (1992b), the similarity between the two samples becomes 
even more apparent (Fks = 16%). 


Fig. 5.— cTg = FWHM[oiii]/2.35 plotted against cr* for the Nelson & Whit¬ 
tle (1996) sample. Open symbols are Type 2 objects; filled symbols are Type 1 
objects. The solid line denotes erg = cr*. The extreme outliers are the powerful 
linear radio sources discussed by Whittle (1992b). 

We can also ask whether Ad is correlated with Eddington 
ratio in Type 1 objects using the Nelson & Whittle (1996) sam¬ 
ple. In Figure 3b we have overplotted the Type 1 objects from 
the Nelson & Whittle sample with measured d*, L[oiii], and 
FWHM[oiii]. Seven of the velocity dispersions (NGC 3227, 
3516, 4051, 4151, 4593, 5548, and Mrk 590) were updated 
from the more recent work of Nelson et al. (2004). We cor¬ 
rected the [O III] luminosity to our assumed cosmology. The 
correlation is not significant for the Type 1 objects (p = 0.27, 
Fnuii = 0.21), possibly because of the small number of objects 
considered, but they seem to trace roughly the trend defined by 
the Type 2 objects. 

5.2. Aperture Ejfects 

Because the SDSS aperture includes the central ~ 5 kpc of 
the galaxy a serious concern in the use of SDSS velocity disper¬ 
sions is that rather than measuring the bulge velocity dispersion 
at a well-defined radius, a distance-dependent fraction of disk 
light may be included in the measurements. Both nuclear gas 
and stellar velocity dispersions may vary with radius. There¬ 
fore, we must examine the extent to which non-nuclear light 
contamination matters. Below we examine both the luminosity 
and kinematic distributions of stars and gas within the SDSS 
aperture to argue that our measurements are meaningful despite 
the large aperture. Rotational broadening is further discussed 
in the Appendix. 

In the case of the gas, we expect the vast majority of the 
emission to come from the NLR at sub-kpc scales. After all, 
our sample was selected to show unambiguous AGN-like line 
ratios. Star-forming regions are the strongest contaminant, * but 
if emission from star formation were important, then the spec¬ 
tra would look like those of H II regions or transition objects 


*We note that of the three tracers considered here, [O ii] is most susceptible to contamination by star formation. Only in low-metallicity H ii regions, which are 
unlikely in the massive hosts of AGNs, is [O iii] a major coolant, while [S ii] is strong in H ii regions only in the presence of shocks (e.g.. Ho et al. 1997c; Kewley et 
al. 2001). 
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rather than bona fide AGNs. We can use the relation between 
NLR radius (Rnlr) and AGN luminosity to estimate the typical 
radius dominating the gas emission. Schmitt et al. (2003), using 
a sample of 60 Seyferts of both types imaged with the Hubble 
Space Telescope, find log Rnlr = (0.33 ± 0.04) x log L[oui] 
- (10.78 ± 1.80), with Rnlr in pc and L[oui] in erg s“^. The 
resulting range in physical size for our objects is ~ 0.25- 1.2 
kpc. At z = 0.1, the corresponding angular size ranges from 
0714 to 0."67, well within the SDSS aperture. Thus, we expect 
the radial variation in (jg to be negligible. 

We now turn to the luminosity distribution of starlight within 
the aperture. Unlike the gas, the stars are distributed through¬ 
out the aperture, and the velocity dispersion we measure may 
differ substantially from cr* at Rg. We estimate the sizes of the 
bulges in our sample, and then infer the probable contamination 
to (T* from disk light. For the median stellar velocity dispersion 
of ~ 150 km s“' in our sample, the fundamental plane relation 
from Bernardi et al. (2003b) predicts « 4.5 kpc. Alterna¬ 
tively, for a median galaxy mass of ^ 5 X 10^° Mq, Kauffmann 
et al. (2003b) find a median Re ~ 3.3 kpc. If the bulges are 
characterized by Re ~ 4 kpc then the majority of the starlight in 
the ~ 5.4 kpc aperture originates in the bulge, and we are not 
probing stars much beyond an effective radius. 

Even if the majority of starlight comes from the bulge, we 
still must consider the radial dependence of ct* with radius. 
Jprgensen, Franx, & Kjaegaard (1995) use 51 E/SO galaxies 
available in the literature to derive an empirical aperture cor¬ 
rection. For E/SO galaxies, the integrated ct* decreases outward 
as a function of radius, falling as a power law with slope -0.04, 
relative to an aperture with a radius of Re/8. On the other hand, 
Pizzella et al. (2004) find that the radial distribution in velocity 
dispersion is flat or slightly increasing outward in a sample of 
spiral (Sb-Sc) galaxies (see their Fig. 5). In early-type galax¬ 
ies a large aperture will tend to underestimate cr* for the bulge, 
while for late-type galaxies there may be a slight tendency in the 
opposite sense. Since our sample consists of early-type galax¬ 
ies, on average the aperture variation most likely causes us to 
underestimate cr*. Bernardi et al. (2003a) find a mean aperture 
correction of ~ 7% for a sample of early-type SDSS galaxies 
based on the empirical function from Jprgensen et al. (1995). 
In any case, the radial distribution of ct* is quite flat. As a 
check, we can look for possible trends in (cTg/cr*) as a function 
of redshift. Since galaxy contamination, on average, should in¬ 
crease as a function of redshift, we might expect (crg/cr*) to 
increase at higher redshift, if we assume (jg = ct* intrinsically. 
We show ((jg/cr*) in three redshift bins (z < 0.1,0.1 < z < 0.2, 
and z > 0.2) in Table 1. We see no significant trend. 

6. CONCLUDING REMARKS 

Using a large sample of Type 2 AGNs selected from SDSS, 
we have systematically investigated the relation between stellar 
velocity dispersions (cr*) and the kinematics of the emission¬ 
line gas ((Jg) as traced by optical forbidden lines ([OII] A3727, 
[O III] A5007, and [S II] AA6716,6731). We show that the 
widths of the low-ionization lines are gravitationally domi¬ 
nated, and may be used as a proxy for cr*. We find that if the 
blue asymmetric wing of the [O III] line is removed, the width 
of the core component can also be used to estimate ct*. The 
ionized gas in the central regions of active galaxies is evidently 
in approximate pressure support with the gravitational potential 
of the bulge stars of the host galaxy, consistent with the situa¬ 
tion in inactive early-type galaxies (e.g., Bertola et al. 1995; 
Cinzano et al. 1999). The scatter in each of the (jg-cr* relations 


is large, however, and we advocate their use only in statistical 
studies. 

While not the primary goal of this study, our results offer 
some insights into the interaction between the AGN and the gas 
of the circumnuclear region. The nearly ubiquitous asymmetric 
blue wing of the [O III] line has long been recognized as indica¬ 
tive of radial motions in the NLR, presumably associated with 
an outflow. In this study, we have shown that this outflowing 
component is principally responsible for imparting supervirial 
motions to the gas. That lines with lower critical density gen¬ 
erally lack this component suggests that it probably originates 
from a more compact region closer to the center. But what is 
the main physical driver of this component? We have shown 
that the most important variable appears to be the Eddington 
ratio. This finding may offer a useful outer boundary condition 
for constraining disk-wind models of the BLR (e.g., Murray & 
Chiang 1997; Proga, Stone, (& Drew 1998; Proga, Stone, (& 
Kallman 2000). 
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APPENDIX 

ROTATIONAL BROADENING 

Rotation must play a role in the nuclear kinematics. In particular, it is not unreasonable to expect some fraction of the nuclear gas 
to reside in a disk-like component. Whittle (1992b) concludes that rotation contributes at least partially to the [O III] line width in 
order to reproduce the observed relations between rotational velocity, inclination, and line width. Vega Beltran et al. (2001) find that 
the ionized gas in a sample of 20 late-type disk galaxies is dominated by rotation. Without rotation velocities or inclination angles for 
our sample, it is difficult to directly examine the rotational contribution. Only the double-peaked emission lines provide evidence that 
rotation can dominate the kinematics, though they are only readily apparent in ~ 1% of the cases. However, we worry that rotational 
broadening may artificially increase our observed line widths. The same concern applies to the stellar velocity dispersions, since the 
SDSS aperture is sufficiently large that some disk light must contaminate the bulge cr*. 

We can estimate the contribution of rotational broadening to the observed scatter from both stellar and gaseous components by 
constructing a toy galaxy model with the average properties of our sample. In particular, we need a model for both the luminosity 
and radial velocity distributions in the disk. The rotational width is then simply 

= J(v-vf I dv/ J I dv, 

where I and v are both functions of radius, and v is the luminosity-weighted mean velocity. Eor the intensity, we use an exponential 
distribution, I = Iq exp (-R/Rj), with an exponential scale length of Rj « 3 kpc (Neistein et al. 1999). Since the gas emission is 
dominated by the NLR, the gas emission radii range from 0.25-1.2 kpc, while the stellar disk extends to a projected aperture radius 
of 2.7 kpc at z = 0.1. 

The maximum circular rotational velocity, Vc, may be estimated using the Tully-Eisher relation. In order to account for the possible 
contamination of the host luminosity from the AGN, we use two different estimates of the absolute magnitude of the underlying host 
galaxy. The first uses the measured, K-corrected values from Brinchmann et al. (2004a), corrected to our assumed cosmology. 
This yields an average = -21.98 mag. One may worry, however, that AGN contamination is present even in Type 2 objects. We 
therefore indirectly estimate in the following way. Erom the relation between concentration index and (M/L)^ shown in Eigure 13 
of Kauffmann et al. (2003a) we can estimate the probable range in (M/L)^. Eor our estimated mean concentration index of ~ 2.65, 
we adopt a range in logarithmic mass-to-light ratios of (M/L)^ = -0.1 to 0.2. If we assume that the typical host galaxy has a mass 
of ~ lO'® ^ Mq, we find a range in absolute magnitude® of -21.74 < M^< -22.49. We are encouraged that the two methods are 
consistent. In what follows, we adopt = -22 mag. Giovanelli et al. (1997) give the following form for the /-band Tully-Eisher 
relation: Mj = -21 -7.52x, where x = log (vJ2) - 2.5. We estimate Mj using the color correction for SO galaxies from Eukugita, 
Shimasaku, & Ichikawa (1995), z-Ic = 0.21 mag. Our resulting average Vc is 230 km s“'. While we argue that our sample is 
dominated by early-type galaxies, we note that recent studies of the Tully-Eisher relation in SO galaxies have found remarkably small 
offsets from the late-type population (Neistein et al. 1999; Hinz, Rieke, & Caldwell 2003). 

As shown by Sofue et al. (1999), the rotation curves of active galaxies are comparable to those of quiescent galaxies, and so we 
use their measured rotation curves of barred and unbarred Sb and Sc galaxies as templates for our model. As a worst-case scenario, 
we only consider edge-on disks with Vc = 230 km s“'. The radius at which the rotation curve flattens is allowed to vary between 0.1 
and 1 kpc, but the rising portion of the rotation curve need not be monotonic; in one of our models the velocity peaks at ^ 0.3 kpc, 
and then dips again before flattening at ~ 2.5 kpc. In the case of the gas, assuming emission radii between 0.25 and 1.2 kpc, Crot 
ranges from 8 to 74 km s“'. The measured a is simply the quadrature sum of the true CTg and arot, so a measured « 150 km s“' 
may be overestimated by up to 15%. To estimate arot in stars, we assume the same family of rotation curves but fix the emission 
radius at 2.7 kpc, for a maximum error of ~ 12% in ct*. If we allow each component to independently assume its maximum error, 
and examine the case where Cg = cr* = 150 km s“', we infer that the observed CTg/cr* may vary from 0.9 to 1.1. We therefore conclude 
that rotational broadening contributes to, but cannot fully account for, the observed scatter. 


®To convert luminosity to absolute magnitude we assume, following Kauffmann et al. (2003a), that the z absolute magnitude of the Sun is q = 4.51 mag at 
z=0.1. 
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